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Abstract 

This paper studies a generalization of the standard continuous-time consensus protocol, obtained by replacing the Laplacian 
matrix of the communication graph with the so-called deformed Laplacian. The deformed Laplacian is a second-degree matrix 
polynomial in the real variable s which reduces to the standard Laplacian for s equal to unity. The stability properties of the 
ensuing deformed consensus protocol are studied in terms of parameter s for some special families of undirected and directed 
graphs, and for arbitrary graph topologies by leveraging the spectral theory of quadratic eigenvalue problems. Examples and 
simulation results are provided to illustrate our theoretical findings. 
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1 Introduction 

In the last decade we have witnessed a spurt of in- 
terest in multi-agent systems research in the control, 
telecommunication and robotics communities (Bullo 
et ai, 2009; Mesbahi and Egerstedt, 2010; Zampieri, 
2008; Kumar et ai, 2008). Distributed control and 
consensus problems (Olfati-Saber et ai, 2007; Ren et 
al., 2007), have had a large share in this research ac- 
tivity. Consensus theory originated from the work of 
Tsitsiklis (Tsitsiklis, 1984), Jadbabaie et al. (Jadbabaie 
et ai, 2003) and Olfati-Saber et al. (Olfati-Saber and 
Murray, 2004), in which the consensus problem was for- 
mulated for the first time in system-theoretical terms. 
A very rich literature emanated from these seminal con- 
tributions in recent years. In particular, numerous ex- 
tensions to tlie prototypal consensus protocol in (Olfati- 
Saber and Murray, 2004) have been proposed both 
in the continuous- and discrete-time domain: among 
them, we limit ourselves to mention here the cases of 
time-varying network topology (Moreau, 2005; Ren 
and Beard, 2005), of networks with delayed (Olfati- 
Saber and Murray, 2004) or quantized/noisy commu- 
nication and link failure (Frasca et ai, 2009; Kar and 
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Moura, 2009), of random networks (Porfiri and Stil- 
wcU, 2007; Tahbaz-Salehi and Jadbabaie, 2008; Fagnani 
and Zampieri, 2008), of networks with antagonistic 
interactions (Altafini, 2012; Altafini, 2013), of dis- 
tributed average tracking (Spanos et ai, 2005; Yang 
et ai, 2008; Chen et ai, 2012), of finite-time con- 
sensus (Wang and Xiao, 2010; Kibangou, 2012), of 
logical (Fagiolini et al., 2008) and nonlinear agree- 
ment (Bauso et ai, 2006; Cortes, 2008), and of consensus 
over finite fields (Pasqualetti et ai, 2013, submitted). 
This paper follows this active line of research and pro- 
poses an original extension to the basic continuous-time 
consensus protocol in (Olfati-Saber and Murray, 2004), 
that exhibits a rich variety of behaviors and whose flex- 
ibility makes it ideal for a broad range of mobile robotic 
applications (e.g., for clustering, or for containment and 
formation control). The new protocol, termed deformed 
consensus protocol, relies on the so-called deformed 
Laplacian matrix, a second-degree matrix polynomial in 
the real variable s, which extends the standard Lapla- 
cian matrix and reduces to it for s equal to unity: the 
deformed Laplacian is indeed an instance of a more 
general theory of deformed differential operators devel- 
oped in mathematical physics in the last three decades, 
cf. (Hislop and Sigal, 1996, Ch. 18). Parameter s has 
a dramatic effect on the stability properties of the de- 
formed consensus protocol, and it can be potentially 
used by a supervisor to dynamically modify the behav- 
ior of the network and trigger different desired agents' 



responses according to time-varying external events. 
The stability properties of the proposed protocol are 
studied in terms of parameter s for some special families 
of undirected and directed graphs for which the eigenval- 
ues and eigenvectors of the negated deformed Laplacian 
can be computed in closed form. In the case of directed 
graphs, it is shown that differently from the standard 
consensus algorithm, for some values of s the states of 
the deformed consensus protocol may also experience 
stable steady-state oscillations. Our analysis is extended 
to arbitrary graph topologies by exploiting the spectral 
theory of quadratic eigenvalue problems (Tisseur and 
Meerbergen, 2001). The discrete-time version of our 
consensus protocol, that involves the so-called deformed 
Perron matrix, is also briefly discussed. 
Beside the aforementioned promising applications, we 
believe that the study of the proposed protocol is of 
value for shedding new light on known results (Olfati- 
Saber and Murray, 2004; Ren and Beard, 2005), and for 
gaining a more general perspective on consensus algo- 
rithms. 

A preliminary version of this paper appeared in (Morbidi, 
2012), compared to which we present here several new 
theoretical results as well as more extensive numerical 
simulations. 

The rest of the article is organized as follows. In Sect. 2 
we review some relevant notions of algebraic graph the- 
ory. The main theoretical results of the paper are pre- 
sented in Sect. 3. In Sect. 4, two possible extensions of 
our results are discussed. Finally, in Sect. 5, the theory 
is illustrated via numerical simulations, and in Sect. 6 
the main contributions of the paper are summarized 
and possible future research directions arc outlined. 

2 Preliminaries 

In this section, we briefly recall some basic notions of 
algebraic graph theory that will be used through the 
paper. Let G = {V, E) be an undirected graplQ where 
V = {1, . . . , n} is the set of vertices, and E is the set of 
edges (Godsil and Royle, 2001). 

Definition 1 (Adjacency matrix A) The adjacency 
matrix A = [oy] oj graph Q is annxn matrix defined as, 

'i if{i,j}eE, 

otherwise. 



Definition 2 (Laplacian matrix L) The Laplacian 
matrix of graph Q is an n x n matrix defined as, 

L = D-A, 



where D = diag{Al.) is the degree matrix^ and 1 = 1„ 
is a column vector of n ones. o 

Note that the Laplacian L is a symmetric positive 
semidefinitc matrix (Mohar, 1991). 

Property 1 (Spectral properties of L) LetAi(L) < 
•^2(L) < ... < A„(L) he the ordered eigenvalues of the 
Laplacian L. Then, we have that (de Ahreu, 2007): 

(1) Ai(L) = with corresponding eigenvector \. The al- 
gebraic multiplicity o/Ai(L) is equal to the number 
of connected components in Q . 

(2) A2(L) > if and only if the graph Q is connected. 
A2 (L) is called the algebraic connectivity or Fiedler 
value of the graph Q . o 

Definition 3 (Bipartite graph) A graph Q is called 
bipartite if its vertex set V can be divided into two disjoint 
sets Vi and V2, such that every edge connects a vertex 
in Vi to one in V2. Equivalently, we have that a graph 
is bipartite if and only if it does not contain cycles of 
odd length. o 

Definition 4 (Signless Laplacian matrix Q) The 

signless Laplacian matrix of graph Q is defined as 
(Cvetkovtc etal., 2007) 

Q = D + A. ^ 

Note that as L, the signless Laplacian Q is a symmetric 
positive semidefinitc matrix (but it is not necessarily 
singular (Brouwer and Haemers, 2012)). Indeed, Q = 
RR^ where R is the vertex-edge incidence matrijff] of Q . 

Property 2 (Spectral properties of Q) The sign- 
less Laplacian Q has the following spectral properties: 

(1) //Ai(L) < A2(L) < ... < A„(L) and Ai(Q) > 
-^2(Q) > ••■ > -^n(Q) fl'^e the ordered eigenval- 
ues of the Laplacian and signless Laplacian, respec- 
tively, then we have that (Cvetkovic et al., 2007): 

A„(L) < Ai(Q). 

Moreover, if n > 2 we have that 

A2(L) < A2(Q) + 2, 

with equality if and only if Q is the complete 
graph Kn (Cvetkovic and Simic, 2009, Th. 3.5). 



^ All graphs in this paper are finite, and with no self-loops 
and multiple edges. 



^ diag(b) is a diagonal matrix with the elements of the vector 
b G H" put on its main diagonal. 

^ The vertex-edge incidence matrix of a graph Q is the 0-1 
matrix R = [vik], with rows indexed by the vertices and 
column indexed by the edges, where rik = 1 when vertex i 
is an endpoint of edge k. 



(2) IfQ is a connected graph withn vertices andm edges, 
then 

, . „ N 2 771 

Ai Q) < -+n-2, 

n — 1 

with equality if and only iJQ is the star graph Ki ^ „_ i 

or the complete graph A'„ (Cvetkovic and Simic, 

2010, Th. 1). 

(3) We have that 

trace(L) = tracc(Q) = 2\E\, 

where \E\ denotes the cardinality of the edge 
set E (Brouwer and Haemers, 2012). 

(4) A graph Q is regular (i.e., each vertex of Q has 
the same degree) if and only if its signless Lapla- 
cian has an eigenvector whose components are all 
ones (Cvetkovic and Simic, 2009, Prop. 2.1). 

(5) IfQ is a regular graph of degree k (i.e., each vertex 
of Q has the same degree k < n — 1), then, 

Pl(A) = (-1)>q(2^-A), 

where pl(A) denotes the characteristic polyno- 
mial of the Laplacian L. If Q is a bipartite graph, 
then (Cvetkovic ei aX., 2007, Prop. 2.3), 

Pl(A) = Pq(A). 

(6) The least eigenvalue o/ Q of a connected graph is 
equal to if and only if the graph is bipartite. In this 
case, is a simple eigenvalue (Cvetkovicct al., 2007, 
Prop. 2.1). 

(7) In any graph, the multiplicity of the eigenvalue 
of Q is equal to the number of bipartite compo- 
nents of the graph Q (Brouwer and Haemers, 2012, 
Prop. 1.3.9). 

(8) Let Q be a regular bipartite graph of degree k . Then 
the spectrum of Q is symmetric with respect to the 
point K (Cvetkovic and Simic, 2009, Prop. 2.2). o 

Let T> = (y, E) be a directed graph (or digraph, for short) 
where T^ = {1, . . . , n} is the set of vertices and E C Vx V 
is the set of edges. In the case of directed graphs, we can 
define the adjacency and degree matrix, as, 

1 if {j, i) e E, 

otherwise, 

and D = diag((ii,i(l), . . . ,din(n)), where din{i) denotes 
the in-degree of vertex i with i G {1, . . . ,n} (i.e., the 
number of directed edges pointing at vertex i). With 
these definitions in hand, the in-degree Laplacian L(2?) 
and in-degree signless Laplacian Q(I') of 2?, can be de- 
fined as in the undirected cas(0- Note that L(I?) and 



Q{T>) are nonsymmetric matrices, and that all the eigen- 
values of L(2?) have non- negative real parts (this can 
be easily proved using Gersgorin's disk theorem (Olfati- 
Saber et ai, 2007)). 

The following definitions will be used in Sect. 4.2. 

Definition 5 (Bipartite digraph) A digraph T> ~ 
(y, E) is called bipartite if its vertex set V can be 
divided into two disjoint sets Vi and V2, such that 
Er\{ViX Vi) = 'and E C^ {V2 x V2) = %, where 
denotes the empty set. O 

Definition 6 (Rooted out-branching) A digraph 
T> = (V, E^ is a rooted out-branching if (Mesbahi and 
Egerstedt, 2010): 

(1) It does not contain a directed cycle; 

(2) It has a vertex ur (root) such that for every other 
vertex V £ V, there is a directed path from vy{ tov.o 

Definition 7 (Strongly connected digraph) A di- 
graph is strongly connected if, between every pair of 
distinct vertices, there is a directed path. o 

Definition 8 (Weakly connected digraph) A di- 
graph is weakly connected if its disoriented version (i.e. 
the graph obtained by replacing all its directed edges with 
undirected ones), is connected. o 

Definition 9 (Balanced digraph) A digraph is called 
balanced if, for every vertex, the in-degree and out-degree 
are equal, i.e., din{i) = rfout(*); for all i (z {1, . . . ,n}. o 

3 Deformed consensus protocol 

3. 1 Problem formulation 

It is well-known (Olfati-Saber and Murray, 2004), that 
if the undirected communication graph Q is connected, 
each component of the state vector x = [xi, . . . , a;„]^ € 
R" of the linear time-invariant system. 



x(i) = -Lx(t), 



(1) 



* "Out-degree" versions of L(0) and Q{T>) can be similarly 
introduced, but they will not be considered in this paper. 



asymptotically converges to the average of the initial 
states xi(0), . . . ,a;„(0), 

I " 1 

lim Xi{t) = - y^ Xi{0) = - xjl, 
t->-oo n ^ — ' n 

i — l 



where Xq = [xi(0), . . . , a;„(0)]"^, i.e., average consensus 
is achieved. The converge rate of the consensus proto- 
col (1) is dictated by the algebraic connectivity A2(L). 
Let us now consider the following generalization of the 
Laplacian L. 
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Fig. 1. Illustrative example: The 6 agents rendezvous at the 
origin while avoiding the two obstacles (gray rectangles): this 
is made possible by switching s from —1 to (the initial 
position of the vehicles is marked with a circle and the final 
position with a diamond) . 

Definition 10 (Deformed Laplacian A(s)) The de- 
formed Laplacian of the graph Q is an n x n matrix 
defined as, 

A(s) = (D-I„)s2 _ As + I„, 

where I„ is the n x n identity matrix, and 
parameter. 



s is a real 
o 



Note that A(s) is a symmetric matrix (but not positive 
semidcfinitc as L, in general), and that: 

A(1)=L, A(-1) = Q. 

Since A(0) = I„, the deformed Laplacian is a comonic 
polynomial matrix (Gohberg et al., 2009, Sect. 7.2). 

Inspired by (1), we will study the stability properties of 
the following linear system. 



x(t) = -A(s)x(t), 



(2) 



in terms of the real parameter s, assuming that the graph 
Q is connected. We will refer to (2), as the deformed 
consensus protocol. 

Remark 1 Note that parameter s in the deformed 
Laplacian A(s) can be regarded as an "exogenous input" 
and it can be exploited to dynamically modify the behav- 
ior of system (2). This may be useful when the vertices 
of the graph are mobile robots and a human supervisor 
is interested in changing the collective behavior of the 
team over time, e.g., by switching from a marginally- to 
an asymptotically-stable equilibrium point of system (2), 
or between two marginally- stable equilibria. The former 
case is illustrated in the example in Fig. 1, where the 
communication graph is the path graph Pg-' in order to 



make the 6 single-integrator agents rendezvous at the 
origin while avoiding the two gray obstacles, the super- 
visor can set s ~ —1 in the first half of the simulation 
and s ~ in the second half (cf. Prop. 1 in Sect. 3.2). o 

Note that since A(l) = L, we will always achieve av- 
erage consensus for s ~ 1. Moreover, since A(s) is 
real symmetric, all the eigenvalues of A(s) (which are 
nonlinear functions of s) are real, and the deformed 
Laplacian admits the spectral decomposition A(s) = 
U(s)A(s)U^(s), where U(s) = [ui(s) U2(s) ... u„(s)] 



R" 



is the matrix consisting of normalized 



and mutually orthogonal eigenvectors of A(s) and 
A(s) == diag(Ai(A(.s)),...,A„(A(s))). The solution 
of (2), can thus be written as, 

x(t) =f^e-^'(^(^))*(ur(s)xo)u,(s). (3) 



In Sect. 3.2, we will focus on some special families of 
undirected graphs for which the eigenvalues and eigen- 
vectors of A(s) can be computed in closed form, and 
thus the stability properties of system (2) can be easily 
deduced from (3). In Sect. 3.3, we will address, instead, 
the more challenging case of undirected graphs with ar- 
bitrary topology. Finally, some extensions (the discrete- 
time case, and the case of directed communication net- 
works), will be discussed in Sect. 4. 

3.2 Stability conditions for special families of graphs 

This section presents a sequence of nine propositions 
which provide stability conditions for system (2), in the 
case of path, cycle, full m-ary tree, wheel, 77i-cube (or 
hypercube), Petersen, complete, complete bipartite and 
star graphs (see Fig. 2 and refer to (Godsil and Royle, 
2001) for a precise definition of these graphs). In the 
following, 

k^ [-1, 1,-1, l,...,(-l)«-i,(-l)«]reIR", 

and hereafter we will use Xi{s) as a shorthand for 
X,{-A{s)),ie{l,...,n}. 

In order to prove our first proposition, we need the fol- 
lowing theorem (Golub and van Loan, 1996, Th. 8.5.1). 

Theorem 1 (Sturm sequence property) Consider 
the following n x n symmetric tridiagonal matrix. 



fli bi 




bi 02 62 




h as 






'■ bn-l 




b„-i an 



where bj 7^ 0, Vj e {1, . . . , n - 1}. Let T^*") denote the 
leading r x r principal submatrix ofT. Then, the number 
of negative eigenvalues of T is equal to the number of 
sign changes in the Sturm sequence; 

1, det(T(i)), det(T(2)), ..., det(T(")). 

The result is still valid if zero determinants are encoun- 
tered along the way, as long as we define a "sign change" 
to mean a transition from + or to ~ , or from — or to 
-\-, but not from + or — to 0. I 



Proposition 1 (Path graph Pn) For the path graph 
P„ with n > 2 vertices (see Fig. 2(a)), we have that: 

• For \s\ < 1, system (2) is asymptotically stable. 

• For \s\ > 1, system (2) is unstable. 

• For s = —I, system (2) is marginally stable. In this 
case, it is possible to identify two groups ofn/2 vertices 
(ifn is even), or one group of [n/2j verticefr] and one 
o/[n/2j +1 vertices (ifn is odd). The states associated 
to the vertices in one group asymptotically converge to 
— xjTlk and the states associated to the vertices in the 



other group converge to —— xjk. 



Proof; In this 
matrix. 


case. 


— A(s) is a symmetric tridi 




' -1 


s 




s 


-(s^ + l) s 


-A(,s) = 




S -(s2 + l) S 






s -1_ 



The Sturm sequence o/— A(s) is given by 

1, —1, 1, —1, 1, . . . , —1, 1, s^ — 1, if n is odd, 
1, —1, 1, —1, 1, ...,1, —1, 1 — s^, if n is even. 

Therefore by Theorem 1, if \s\ < 1, s 7^ 0, all the eigen- 
values of — A(s) are strictly negative and system (2) is 
asymptotically stable. On the other hand, the Sturm se- 
quence of ^{s), for all n, is given by 1, 1,..., 1, 1 — s^, 
from which we deduce that for \s\ > 1, A(.s) has a nega- 
tive eigenvalue, and hence system (2) is unstable. Since 
A(0) = In, the system is asymptotically stable for s = 0. 
Finally, for s = —1, system (2) is marginally stable and 
the unit-norm eigenvector associated to the zero eigen- 
value of -A(-l) is -^k. ■ 



^ [-J denotes the floor function, which maps a real number 
to the largest previous integer. 



Proposition 2 (Cycle graph C„) For the cycle graph 
C„ with n > 2 vertices (see Fig. 2(b)), we have that: 

• If n is even: 

■ For all s G ]R\{ — 1, 1}, system (2) is asymptotically 
stable. 

■ For s = —\, system (2) is marginally stable. In this 
case, the states associated to n/2 vertices asymptoti- 
cally converge to — Xq k and the states associated to 
the other n/2 vertices converge to — Xq k. 

• If n is odd, system (2) is asymptotically stable for all 
sgR\{1}. 

Proof; In this case — A(s) is a circulant matrix, 
-A(s) == circ[-(s2 + 1), s, 0, . . . , 0, s], 

i.e., each subsequent row is simply the row above shifted 
one element to the right (and wrapped around, i.e., mod- 
ulo n). The entire matrix is thus determined by the first 
row. It is well-known that circulant matrices are diago- 
nalizable by the Fourier matrix F„, given via. 



1 1 

1 UJ 



1 



1 



UJ 



1 

,,2(«-l) 



1 w"^-^ aj^("^^) • • • w("-i)("^i) 



where u = g^'^^'", j = \/— 1, and F* denotes the con- 
jugate transpose o/F„, hence their eigenvalues can be 
computed in closed form. The eigenvalues of a general 
n X n circulant matrix C ~ circ[ci, C2, . . . , c„], in fact, 
are given by: 



UC) = pc{lo'-'), ze {!,..., n}, 



(4) 



where the polynomial pc{C) = CnC^^^ + . . . + cs^^ + 
C2^+ Ci is called the circulant's representer (Davis, 1994, 
Th. 3.2.2). By applying this result to matrix — A(s), for 
i £ {1, . . . , n} we have that, 



A.(-A(s)) 



2 cos 



2 7r(i-l) 



s - 1. (5) 



Observe now that the coordinates of the vertex of the 
parabola (5) are 



cos 



2Tr{i - 1) 



2Tr{i - 1) 



, i e {1, ..., n}. 



// n is even, then Ai(— A(s)) < 0, Vs G R and 
\/i ^ {1, n/2 + 1}. For i = 1, Ai(-A(s)) < and 
Ai(-A(s)) = only for s = 1. For i = n/2 + I, 



@ @ 
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(e) 



(f) 



(g) 



(h) 
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(i) 



Fig. 2. Families of undirected graphs: (a) Path graph P„; (b) Cycle graph Cn', (c) FuU m-ary tree; (d) Wheel graph Wn] 
(e) m-cube Qm {m = 4); (f) Petersen graph; (g) Complete graph Kn\ (h) Complete bipartite graph Km,n\ (i) Star graph K\^n- 

\^i2 + i{—^{sy) < and A„/2 + i(— A(s)) = only levels of the tree asymptotically converge to -i=XqUi, 



for s = —1. The unit-norm eigenvector associated to 
A„/2 + i(^A(^l)) *■? -7= k. On the other hand, if n is 
odd, then X^{- A{s)) < 0, V.s e II andVi ^ 1. Fori = 1, 
Ai(-A(s)) <0 andXi{-A{s)) = only for s ^ 1. ■ 



Note that 



det(-A(s)) ^ (-l)"(s2 



2s" 



1), 



for the cycle graph C„, and that the 2n roots of 
det(— A(s)) are evenly spaced on the unit circle. 

A full m-ary tree is a rooted tree in which every ver- 
tex other than the leaves has m children (2-ary and 3- 
ary trees are sometimes called binary and ternary trees, 
respectively (Godsil and Royle, 2001)). The depth 6 of 
a vertex is the length of the path from the root to the 
vertex. The set of all vertices at a given depth is called 
a level of the tree: by definition, the root vertex is at 
depth zero. The number of vertices of a full m-avy tree 

is n = X]i=o"^'- 

Proposition 3 (Full in-ary tree) For the full m-ary 
tree, with m >2 (see Fig. 2(c)), we have that: 

• For \s\ < 1, system (2) is asymptotically stable. 

• For \s\ > 1, system (2) is unstable. 

• For s = — 1 , system (2) is marginally stable. In this 



while the states associated to the vertices in the odd 



levels of the tree converge to 



\/n 



Xg Ui , where Ui is 



the unit-norm eigenvector associated to the zero eigen- 
value o/^A(— 1). 

Proof; The stability properties of (2) are determined in 
this case by only one of the eigenvalues of — A(s) (in 
fact, the other n — 1 are negative for all s: note that 
det(— A(s)) = s^ — 1). This eigenvalue is negative for 
\s\ < 1, positive for \s\ > 1 and zero for s e { — 1, +!}• ■ 

Proposition 4 (Wheel graph Wn) Consider a wheel 
graph Wn with n > 3 vertices where vertex 1 is the center 
of the wheel (see Fig. 2(d)), and let fi be the non-unitary 
root of 



yJ{{n~4)s + 2Y +4{n~l) 



1. (6) 



case, the states associated to the vertices in the 



even 



fj. monotonically decreases from 1/2 (for n = 4) to 
(for n = 00 j [see Fig. 3]. We have that: 

• For s > \ or 3 < ^, system (2) is asymptotically stable. 

• For s G (/.J, 1), system (2) is unstable. 

• For s ~ fi, system (2) is marginally stable. In this 
case the state associated to vertex 1 asymptotically 
converges to X(^[a^, a/3, ...,a/3]"^, and the states 
associated to the other n — 1 vertices converge to 
xJ[a/3, /32, . . . , I3Y, where [a, /3, . . . , (3^ , a, /3 e R, 



0.45 



0.35 




0.25 



100 



Fig. 3. Value of the non-unitary root of (6) for n G [4, 100]. 

is the unit-norm eigenvector associated to the zero 
eigenvalue o/ — A(/i). 

Proof.' The eigenvalues of matrix — A(s) are: 



• For s E {—1, tt}; system (2) is marginally sta- 
ble. In this case, the states associated to n/2 ver- 
tices asymptotically converge to -7=XqUi, while the 

states associated to the other n/2 vertices converge 
to — -^XqUi, where Ui is the unit-norm eigenvec- 
tor associated to the zero eigenvalue of — A(— 1) or 

Proof.' In this case, the eigenvalues of —A.{s) are: 

Al(s) = -(("1- l)s2 -1713 + 1), 

his) = -((m- l)s^ + ms + 1), 

A3(s) = ... = A^^2(«) = 

-((m - 1) s^ - (m - 2) s + 1), 



An^ 3(5) = ... = A„(s) 



-((m - 1) s^ + (to - 2) s + 1). 



Ms)^-^s'+s- 



Hs) = -^s^+s- 



y(("- 


-4)s + 2)%4(n- 


-1) 


2 


\/(("- 


-4).s + 2)%4(n- 


-1) 



A,+i(s) = -2s^ + 2cos{i^^^)s-l, z e{2, 



!}• 



A^ote t/iat t/ie coordinates of the vertex of the parabola 
Xi+iis), i e{2,...,ri- 1}, are 



1 /27r(i-l) 
- cos 

2 V n- 1 



1 + sin 

2V V n-1 



therefore Xi+i{s) < 0, Vs € R, Vi £ {2,...,ri- 1}. 
W^e aZso have that A2(s) < 0, Vs e R. Finally, it is 
easy to verify that Ai(s) has always two roots, s — fi and 
s = 1. Ai(s) > for s G (/i, 1) and Ai(s) < /or s > 1 
or s < fi. ■ 

Proposition 5_(m-cube Qm) For the m-cube (or hy- 
percube) grapl\2\ Qm with n = 2™ > 4 vertices (see 
Fig. 2(e)), we have that: 

• For \s\ > 1 or\s\ < — ^, system (2) is asymptotically 
stable. 

• Fors e (- 1, --^) orse {:^, l), system (2) is 
unstable. 

• Fors — ^_^ , average consensus is achieved. The con- 
vergence rate to ^ X(^l„ is slower for s = ^_^_^ than 
for s = 1. 



m-cube or liypercube graphs should not be confused with 
cubic graphs which are 3-regular graphs. The only hypercube 
which is a cubic graph is Qs- 



an I 



Since \^{s), . . ., A„(s) are negative for all s, the stability 
of system (2) is only determined by the eigenvalues Ai(s) 
and\2{s). We have that \i{s) < for s > 1 or s < ^_^ . 

Instead, A2(s) < for s > — ^_^^_ ^ or s < —1. Hence, for 
> 1 or \s\ < ^\-^ , system (2) is asymptotically stable 
d for s € {—1, =fc ^j_.^ }, marginal stability is achieved. 

Finally, note that A.( ^—r) = —^ A(l) = —^ L which 

^' \m— 1/ m — 1 ^ ' m — 1 

implies that \i{^^{^ — ~r)) ^ ^-, A2 (A(l)). //u;e now 
recall the role played by the algebraic connectivity of the 
graph in the consensus protocol, we have that the con- 
vergence rate to — xjl„ is slower for s — — ^ than 
for s = 1. M 

Proposition 6 (Petersen graph) For the Petersen 
graph (see Fig. 2(f)), we have that: 

• For s > 1 or s < 1/2, system (2) is asymptotically 
stable. 

• For s G (1/2, 1), system (2) is unstable. 
For s ~ 1/2, average consensus is achieved. The con 

To 
for s — \ 



vcrgence rate to ^irx^lio is slower for s — 1/2 than 



Proof.' The ten eigenvalues of —A.{s) are: 
Ai(s) == -{2s^-3s + l), 



Hs) 
Xe{s) 



= his) = -i2s^-s + l), 
= Aio(s) = -(2s2 + 2s-f 1). 



We have that A2(s) < and Xq(s) < 0, Vs e R. 
Moreover, Ai(s) < for s > 1 or s < 1/2, and the 
unit-norm eigenvector associated to Ai(l/2) is — == Iiq. 



Graph name 


Asymptotic stability for : 


Marginal stability for : 


Path graph P„,n>2 


\s\<l 


s — —1 (2 groups of vertices) 


Cycle graph Cn, n > 2, n even 


VsG ]R\{-1, 1} 


s — —1 (2 groups of vertices) 


Cycle graph C„, n > 2, n odd 


VsG R\{1} 




Full m-ary tree, m > 2 


ls|<l 


s — —1 (2 groups of vertices) 


Wheel graph W„, n > 3 


s > 1 or s < n 


s — ^ [2 groups of vertices) 


m-cube, Qm, n = 2™ > 4 


\s\ > 1 or \s\ < J_^ 


s G { — 1, — ;;7^} (2 groups of vertices) 
s = ^^^ (average consensus) 


Petersen graph 


s > 1 or s < 1/2 


s = 1/2 (average consensus) 


Complete graph K„, n > 2 


s > 1 or s < -^ 


s = j^ (average consensus) 


Complete bipartite graph 
Km^n, m,n>2 




-^(m-l)(n-l) 

-^ (if m = n, average consensus) 
i/(m-i)(n-i) (jf m^ n, 2 groups of vertices) 


' ' ' y(m-l)(n-l) 


Star graph Ki^n, n > 3 


Is] < 1 


s — —1 (2 groups of vertices) 



Table 1 

Summary of the stability properties of the deformed consensus protocol (2), for some special families of undirected graphs. 

Average consensus is achieved in all cases for s = 1 



Finally, A(l/2) = i A(l), hence the convergence rate 
to -^XqIio is slower for s = 1/2 than for s =1. I 

Proposition 7 (Complete graph Kn) For the com- 
plete graph A'„ with n > 2 vertices (see Fig. 2(g)), we 
have that: 

• For s > 1 or s < —^ , system (2) is asymptotically 
stable. 

• For s G (^ro' ^) ' •''V^^^''^ (^) ''■^ unstable. 

• For s = ;;32' '^^e'''^ffc consensus is achieved. The con- 
vergence rate to ^ Xq 1 is slower for s = ^^^ than 
for s = 1. 

Proof; The eigenvalues o/ — A(s) are: 

Xi{s) = -{{n-2)s^-{n-l)s + l), 



• Fors G (—1, 



, = or s G —^ 



-l)(n-l)' ^' 



.system (2) is unstable. 
• For s G I — 1, ± , ^ , I, system (2) is 

marginally stable. In particular, given the initial con- 
dition ^^ £ Wr+'' : 

■ If m 7^ n: for s = —1, the .states associated 
to the vertices in Vi asymptotically converge to 
^ Xq[1L^, — 1^]"^ and the states associated to the 



m-\-n 



vertices in V2 converge to — „^ , ,^ Xq [1^, 



-1 



T]T 



For s ~ ±- 



--, the states associated to the 



A2(s) 



A„(.s) 



-{{n-2)s^ + s + l). 



We have that A2(s) < 0, Vs G R. Moreover, Ai(s) < 
for s > 1 or s < —z^, and the unit-norm eigenvec- 
tor associated to Ai(— ^2) is -7=1. Finally, note that 
A(— ^2) — "^2^(1)' hence, the convergence rate to 
— xiTl is slower for s = —^ than for s = 1. ■ 

71 ^ •' 71—2 •' 

Proposition 8 (Complete bipartite graph A'™ „) 

For the complete bipartite graph Km,7i = {Vi U V2,E), 
where \Vi\ = m, IV2I = n with m,n > 2 (see Fig. 2(h)), 
we have that: 



V(m-l)(Ti-l)' 

vertices in Vi and V2 asymptotically converge to one 
of the two different values taken by the components 
of the vector (ui u^)xo, where Ui is the unit-norm 
eigenvector associated to the zero eigenvalue of 

'^ y/(7n-l)(7i-l)' 

If 771 = n: for s = ^^^ average consensus is 
achieved, and the convergence rate to ^ x^l2n 
is slower for s = —^ than for s = 1. For 
s £ {~ 1; rrl the states associated to the vertices 

in Vi asymptotically converge to 2;^ Xq [1^, — 1^]^ 
and the states associated to the vertices in V2 con- 
verge to -^x^[ll, -1^]^. 



Proof; In this case, the eigenvalues o/— A(s) are: 



Ai(5) 



n + m — 2 o \/(n — m)^s-^ + Amn 
7: s + -^ '- s - 1, 



• _For |s| > 1 or |s| < 
totically stable. 



^(m-l)(n-l) 



, system (2) is asymp- 



, , , n + m — 2 o \/(n — m)'^s'^ + Amn 

Ms) = ^ s ^ s - 1, 



X,{.s) = ... = A„+i(s) = -{{m - 1) s^ + 1), 

A„+2(s) = . . . = A„+„i(s) = -{{n - 1) s^ + 1). 

Since A3(s), . . ., A„+m(s) are negative for all s, the sta- 
bility of system (2) is only determined by the eigenval- 
ues Ai(s) and A2(s). A systematic study of the roots of 
Ai(s) and A2(s) immediately leads to the result. Note 
that if m = n, A(— ^j) = — ^ ^(1)? 0''^d hence the rate 
of convergence to ^ x^l2n is slower for s — -^^ than 
for s = 1 . I 



From Prop. 8, we can deduce the following result 
(note that the star graph is a complete bipartite graph 
with TO = 1): 

Proposition 9 (Star graph Ki^) For the star graph 
Ki_n with n > 3, where vertex 1 is the center of the star 
(see Fig. 2(i)), we have that: 

• For \s\ < 1, system (2) is asymptotically stable. 

• For \s\ > 1, system (2) is unstable. 

• For s = —1, system, (2) is marginally stable. In 
this case, the state associated to vertex 1 asymp- 
totically converges to ^^^ Xq [1, — l^]-'" and the 
states associated to the other n vertices converge to 

n+l ^0 [-•-' -"^nJ • ■ 

For the reader's convenience, all the results found in this 
section are summarized in Table 1. 

Remark 2 Note that the path, cycle (with n even), full 
m-ary tree, m-cube, complete bipartite and star graphs 
are all bipartite graphs. Then, in view of Property 2. 6, 
for s = — 1 the state of system (2) asymptotically con- 
verges to (u„u^)xo with these graphs, where u„ is the 
unit-norm eigenvector associated to the zero eigenvalue 
of the signless Laplacian Q. According to (Altafini, 2013, 
Def. 1), we can say that we achieve "bipartite consensus" 
in these cases. o 



3.3 Stability conditions for graphs of arbitrary topology 

In order to extend the analysis of the previous section 
to arbitrary undirected graphs, we briefly review here 
the spectral theory of quadratic eigenvalues problems 
(QEPs) (Tisseur and Meerbergen, 2001, Sect. 3), that 
constitute an important class of nonlinear eigenvalue 
problems. Let 

P(A) = BaA^ + BiA + Bo, 

be an n X n matrix polynomial of degree 2, where B2, 
Bi and Bq are n x n complex matrices (Gohberg et 
at, 2009). In other words, the components of the matrix 
P(A) are quadratic polynomials in the variable A. 



Definition 11 (Spectrum of P(A)) The spectrum of 
P(A), denoted by S(P), is defined as, 

I](P) = { A G C : det(P(A)) = 0}, 

i.e., it is the set of eigenvalues 0/ P(A) (Tisseur and 
Meerbergen, 2001). O 

Definition 12 (Regular P(A)) The matrix P(A) is 
called regular when det(P(A)) is not identically zero for 
all values of X, and nonregular otherwise (Tisseur and 
Meerbergen, 2001). O 

Note that det(P(A)) = det(B2)A2"+ lower-order terms, 
so when B2 is nonsingular, P(A) is regular and has 2n 
finite eigenvalues. When B2 is singular, the degree of 
dct(P(A)) is r < 2ri andP(A) has r finite eigenvalues and 
2n — r infinite eigenvalucfQ. The algebraic multiplicity 
of an eigenvalue Aq is the order of the corresponding zero 
in det(P(A)), while the geometric multiplicity of Aq is 
the dimension of ker(P(Ao)). 

Problem 1 (Quadratic eigenvalue problem, QEP) 

The QEP consists of finding scalars A and nonzero 
vectors z, y, satisfying (Tisseur and Meerbergen, 2001), 



P(A)z = 0, y*P(A) 



0, 



where z, y e C" are respectively the right and left eigen- 
vector corresponding to the eigenvalue A € C, and y* is 
the conjugate transpose o/y. <> 

A QEP has 2n eigenvalues (finite or infinite) with up to 
2n right and 2n left eigenvectors. Note that a regular 
P(A) may possess two distinct eigenvalues having the 
same eigenvector. If a regular P(A) has 2n distinct eigen- 
values, then there exists a set of n linearly independent 
eigenvectors. 

Property 3 (Spectral properties of P(A)) // ma- 
trices B2, Bi, Bq are real symmetric, the eigenvalues 
of P(A) are either real or occur in complex- conjugate 
pairs, and the sets of left and right eigenvectors coin- 
cide (Tisseur and Meerbergen, 2001). o 

By leveraging the previous facts, we deduce the following 
property of the deformed consensus protocol (2): 

Proposition 10 The finite real eigenvalues A of 
the QEP, 



((I„-D)A2+AA-I„)z 



0, 



(7) 



are the values of s for which system (2) is marginally 
stable. Moreover, if A is one of these eigenvalues with 



^ The infinite eigenvalues correspond to the zero eigenvalues 
of the reverse polynomial A^ P(A~^) = Bo A^ -f Bi A -I- B2. 



geometric multiplicity 1, andz = z/||z|| 
unit-norm, eigenvector, we have that: 



the associated 



lim x(t) 

t->-oo 



(zz^ 



Ixq. 



Remark 3 (Computation of the eigenvalues) The 

eigenvalues of the QEP (7) can be easily computed by 
converting it to a standard generalized eigenvalue prob- 
IcnPI of size 2n, by defining the new vector w = Az. In 
terms ofz andw, problem (7) becomes: 



' 


In ' 




z 


= A 


[l„ 


A 




w 





I« 

I„ D 



Matlab's function "polyeig" uses a linearization proce- 
dure similar to the one that we have described, to nu- 
merically solve generic polynomial eigenvalue problems. 
Note that if matrix I„ — D is singular, it is more con- 
venient from a numerical viewpoint to use the second 
instead of the first companion form considered above, 
for the linearization (cf. (Tisseur and Meerbergen, 2001, 
Sect. 3.4)). o 

Tiie following proposition elucidates the connection ex- 
isting between the topology of the communication graph 
Q, and the properties of the QEP (7). 

Proposition 11 If the graph Q has a vertex with degree 
equal to one (i.e. only one edge is incident to that vertex), 
the matrix (I„ — D)A^ + A A — I„ is nonregular and the 
QEP (7) admits at least two infinite eigenvalues. ■ 

For example, with the path graph rank(I„ — D) — n — 2, 
with the star graph rank(I„ — D) = 1 and with the full 
m-aiy tree graph rank(I„ — D) = n — m . in the first two 
cases det((I„ - 'D)\^ + AA - I„) = (-1)"-1(A2 - 1), 
while in the last det((I„ - 'D)\^ + A A - I„) = A^ - 1. 
Hence, in all cases, the QEP (7) admits 2n — 2 infinite 
eigenvalues. 

The next proposition, the last result of this section, 
shows how to determine the s-stability interval of the 
deformed consensus protocol, for graphs with arbitrary 
topology (sec Fig. 4). 

Proposition 12 (Stability interval for arbitrary Q) 

Let q{s) = det((I„ - D) s^ + A s - I„), then: 

• Ifn is even, system (2) is asymptotically stable for all 
s such that q{s) > 0, and unstable for all s such that 
q{s) < 0. 

• If n is odd, system (2) is asymptotically stable for all 
s such that q{s) < 0, and unstable for all s such that 
q{s) > 0. 



This construction is called "linearization" in the litera- 
ture and it is not unique in general (see (Demmel, 1997, 
Sect. 4.5)). 



unstab. 




Fig. 4. Illustration of Prop. 12: For the path graph Pq, 
q{s) = l-s^. 

Proof.- From Property 3 and Property 1.1, it follows that: 

i r-2i-l 

9(.s) = (s-i)nG^-a)(s-a) n (^-^^O' 



fc=l 



J = l 



where r < 2n is the degree of q{s), (^fe, ^^) ai^e the 
i pairs of complex-conjugate roots of q{s) and rjj the 
r—2(. — 1 non-unitary real roots ofq{s) (counted with their 

multiplicity). Since Y[k=ii^ ~ Ck){s — Ck) ~ Y[k=i *^ ~ 
23fi[Cfc] s -|- ICfeP > for all s, where 5R[ • ] denotes the real 
part of a complex number, from Prop. 10 the statement 
follows. M 

The following example illustrates the rich variety of be- 
haviors exhibited by the deformed consensus protocol on 
four "generic" (nonbipartite) graphs with five vertices. 

Example 1 Consider the four graphs reported in Fig. 5. 
By leveraging Prop. 10 and Prop. 12, we have that: 

• With the graph in Fig. 5(a), system (2) is asymptoti- 
cally stable y s S ]R\{1}. 

• With the graph in Fig. 5(b), system (2) is asymptoti- 
cally stable for s < 0.7022 or s > l.Fors = 0.7022, the 
system is marginally stable and three groups of vertices 
can be identified: {!}, {2,5}, {3,4} (different shapes 
are used in Fig. 5(b) to indicate these groups). 

• With the graph in Fig. 5(c), system (2) is asymptoti- 
cally stable for s < 0.4396 or s > 1. For s = 0.4396, 
the .system is marginally stable and three groups of ver- 
tices can be identified: {1}, {2,5}, {3,4}. 

• With the graph in Fig. 5(d), system (2) is asymptoti- 
cally stable for s < 0.3804 or s > 1. For s = 0.3804, 
the system is marginally stable and two groups of ver- 
tices can be identified: {1,3}, {2,4,5}. 

From Figs. 5(b)- (d), we notice that vertices in the same 
group tend to have the same edge degree, and that an 
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(a) 



(b) 




(c) 




(d) 



Fig. 5. Example 1: In (b)-(d), difTerent shapes are used to identify distinct groups of vertices: the states associated to the 
vertices in these groups asymptotically converge to the same value when system (2) is marginally stable. 



increase in the algebraic connectivity of the graph leads 
to a reduction of the s-stability interval of the deformed 
consensus protocol. In future works, we are going yo 
study whether the notion of nontrivial equitable parti- 
tion (Mesbahi and Egerstedt, 2010, p. 34), may help to 
shed some light into the peculiar grouping behavior exhib- 
ited by the vertices of the four graphs considered in this 
example. 

4 Extensions 

In this section, two extensions to the theory presented 
in Sect. 3 are discussed: we first consider a discrete-time 
version of the deformed consensus protocol (2) , and then 
deal with directed communication graphs. 

4.1 Discrete-time deformed consensus protocol 

Following (Olfati-Saber et al., 2007, Sect. IIC), we can 
introduce the following discrete-time version of proto- 
col (2), 

x(fc + l) = P(s)x(fc), fcG{0, 1, 2,...}, (8) 



where 

P(s) =I„-6A(s), 

is the deformed Perron matrix and < e < 



(E 



j#* "-y 



,) is the 



is the step-size, where dmax = i 
maximum degree of Q. It is easy to verify that the 
continuous-time and discrete-time deformed consensus 
protocols share the same s-stability intervals, and thus 
the analysis of Sect. 3, mutatis mutandis, is still valid 
for system (8). 

4.2 Deformed consensus protocol for directed graphs 

In this section we assume that the communication graph 
is directed and contains a rooted out-branching, and by 
mimicking Sect. 3 we will study the stability properties 
of the following linear system. 



i(t) = -A(I?(s))x(0, 



(9) 



in terms of the real parameter s, where the symbol 
"I?(s)" indicates that the deformed Laplacian is now 



relative to a directed communication topology. Similarly 
to Sect. 3.1, we have here that: 

A(P(1)) - Lp), A(P(-1)) = Q(P). 

It is well known (Mesbahi and Egerstedt, 2010), that if 
the digraph V contains a rooted out-branching, the state 
trajectory of system. 



satisfies. 



x(i) = -L(2?)x(t), 



lim x(t) = (uivf)xo. 



(10) 



where Ui and vi, are, respectively, the right and left 
eigenvectors associated with the zero eigenvalue of L(I?) , 
normalized such that u^Vi = 1. Moreover, we have 
that (10) reaches average consensus for every initial state 
if and only if V is weakly connected and balanced. 
In what follows, we will analyze the stability properties 
of system (9) for two special families of directed graphs, 
and briefly explore the case of digraphs of arbitrary 
topology with the help of few significative examples. 

Proposition 13 (Directed path) For the directed 
path graph with n > 2 vertices (see Fig. 6(a)), we have 
that: 

• For \s\ < I, system (9) is asymptotically stable. 

• For \s\ > 1, system (9) is unstable. 

• For s G { — 1, +1}, system (9) is marginally stable. 




® — >(D — KD 

(a) 



Fig. 6. Two families of directed graphs: (a) Directed path; 
(b) Directed cycle. Note that the former is weakly connected 
and not balanced, and the latter is strongly connected and 
balanced. 
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• For s ~ —1, it is possible to identify two groups of 
n/2 vertices (ifn is even), or one group of ln/2\ ver- 
tices and one of \n/2\ + 1 vertices (ifn is odd). The 
states associated to the vertices in one group reach 
an agreement on a;i(0), and the states associated to 
the vertices in the other group agree on — a;i(0). 

• For s = 1, the consensus value is a;i(0). 

Proof; In this case, — A(I?(s)) is a lower-triangular 
matrix, 



-A(P(,s)) 



s'-l 
s -1 



s -1 

s -1 



whose eigenvalues are simply, 

Xi{s) = s^-l, A2(s)= ... =A„(.s) = -l, 

from which the statement immediately follows. M 

Protocol (9) exhibits a larger pool of behaviors with the 
directed cycle than with the directed path, as illustrated 
in the next proposition. 

Proposition 14 (Directed cycle D„) For the di- 
rected cycle graph Dn with n > 2 vertices (see Fig. 6(b)), 
we have that: 

• For s — 1 and for all n, average consensus is achieved. 

• If n is even: 

■ For \s\ < 1, .system (9) is asymptotically stable. 

■ For \s\ > 1, system (9) is unstable. 

■ For s = — 1, system (9) is marginally stable. In this 
case, the states associated to n/2 vertices asymptoti- 
cally converge to — Xq Ik and the states associated to 
the other n/2 vertices converge to Xq k. 

• If n is odd: 

■ System (9) is asymptotically stable for s G {'&{n), 1), 
where, 



d{n) 



1 



n(n - 2) + 1 

cos I TT 



For s < ??(n) or s > 1, system (9) is unstable. 
For s = i?(n), system (9) is marginally stable. 
At steady-state, we have that the i-th state of sys- 
tem (9), i £ {1, . . . , n}, obeys, 

x^{t) = Asin(27r/(n)f + 0j(n) + 0o), 

where A and 0o are positive constants, the frequency 

/(n) = -tan . , 



and for i e { 1 , . . . , n} , the phase 
271 {i- 1) 



n tani ^^ — '-^ tt 



4),{n) 



Proof.' Similarly to Prop. 2, — A(2?(s)) is a circulant 
matrix in this case, 

-A{V{s)) = drc^l, 0, . . . , 0, s], 

and using formula (4) the eigenvalues of —A.{T>{s)) can 
be computed in closed-form, 



A,(,s) = . cos( ^-<-;)<"-^> ) 



- 1 



J s sm 



27r(z-l)(n-l) 



j , i e {!,..., n}, 



(11) 



where j is the imaginary unit. The statement easily fol- 
lows from a systematic study of (11) in terms of param- 
eter s. Note, in particular, that ifn is odd and s ~ '&{n), 
the following pair of purely imaginary eigenvalues ap- 
pears in (11), 



, n(n-2) + 1 
Ar„/2l, rr,/2l+i == ±j tanI TT 



where [•] denotes the ceiling function. These eigenvalues 
are responsible for the periodic solutions of period. 



1 



27r 



f{n) tan(^i(^^^7r)' 

of system (9), (cf. (Strogatz, 1994, P- 134)). ■ 

It is easy to verify by direct calculation, that, 

i?(3) = -2, lim din) = -1, lim f{n) = 0, 

n—>-oo n— ^oo 

and that 

dct(-A(D„(s))) = (-l)«(l-s"), 

whose roots are evenly spaced on the unit circle. Note 
that the directed path and the directed cycle (for n even), 
are bipartite digraphs (cf. Remark 2). 

For easiness of reference, the results found in this section 
are summarized in Table 2. 

Remark 4 Note that differently from Sect. 3, since 
A.(T>{s)) is nonsymmetric it can also admit complex- 
conjugate eigenvalues, and the states of system (9) can 
experience stable steady-state oscillations (as we have 
seen in Prop. 14 for n odd). It is worth pointing out 
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Digraph name 


Asymptotic stability for : 


Marginal stability for : 


Directed patli, n > 2 


1^1 <1 


s — —1 (2 groups of vertices) 
s = 1 (consensus) 


Directed cycle Dn, n > 2, n even 


Is] < 1 


s — —1 (2 groups of vertices) 
s = 1 (average consensus) 


Directed cycle D„, n > 2, n odd 


s e (i9(n), 1) 


s = i9(n) (stable oscillations) 
s = 1 (average consensus) 



Table 2 

Summary of the stability properties of the deformed consensus protocol (9), for two families of directed graphs. 







(a) 



(b) 



(c) 



(d) 



Fig. 7. Example 2: Variations on the directed cycle with five vertices. 



here that this is not true, instead, for protocol (10), 
which does not admit periodic solutions (cf. Prop. 3.10 
in (Mesbahi and Egerstedt, 2010)). o 

Since the study of the behavior of system (9) for digraphs 
T) of arbitrary topology is nontrivial, we will focus here 
on few representative examples and try to deduce some 
criteria of general validity. 



Example 2 

in Fig. 7. 



Consider the Jour digraphs reported 



With the digraph in Fig. 1(a), we have that system (9) 
is asymptotically stable for s G (1/ cos(167r/5), 1) 
(cf. Prop. 14). For s = 1/ cos(167r/5), system (9) is 
marginally stable: at steady-state, its states oscillate 
with the same frequency and amplitude, and the phases 
are evenly spaced. For s = 1, average consensus is 
achieved (see Sect. 5.2). 

With the digraph in Fig. 7(b), system (9) is asymptot- 
ically stable for s £ (-1.6889, 1). For s = -1.6889, 
system (9) is marginally stable. At steady-state, its 
states oscillate with the same frequency but have dif- 
ferent amplitudes: the phases are regularly spaced. For 
s = 1, consensus (but not average consensus, since 
edge (1, 3) breaks the balancedness of the original 
five-vertex directed cycle), is achieved (see Sect. 5.2). 
With the digraph in Fig. 7(c), we have the same qual- 
itative behavior as with the digraph in Fig. 7(b). The 
only difference is represented by the stability threshold, 
that is now s = —1.9441 instead of s — —1.6889. 



• With the digraph in Fig. 7(d), system (9) is asymptot- 
ically stable for s £ (—1.3326, 1). For s = —1.3326, 
system (9) is marginally stable and its non-oscillating 
states asymptotically converge to five different val- 
ues. For s — 1, consensus (but, again, not average 
consensus) is achieved. 

Note that we cannot leverage Prop. 12 to determine the 
stability interval of system (9) when its states experi- 
ence stable oscillations (the threshold values for the di- 
graphs in Figs. 7(b) and 7(c), have been determined on a 
trial- and- error basis): however. Prop. 12 appears to pro- 
vide the correct stability intervals in all the other cases 
(e.g., for the directed path in Fig. 6(a), or for the di- 
graph in Fig. 7(d)), and in those cases Prop. 10 seems 
to be valid as well. The analytical determination of these 
stability thresholds for general digraphs is the subject of 
on- going research. The directed graphs in Figs. 7(b) and 
7(c) are particularly significative for showing the non- 
trivial connection existing between the topology of the di- 
graph T) and the threshold values: in fact, here, a single 
edge- orientation change {i.e., (1, 3) versus (3, 1)), yields 
two remarkably different thresholds: s = —1.6889 and 
s = -1.9441. o 

5 Simulation results 

Numerical simulations have been performed in order to 
illustrate the theory presented in Sect. 3 and Sect. 4.2. 
Sect. 5.1 deals with the case of undirected graphs, and 
Sect. 5.2 with the case of directed graphs. 
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t = 1 sec. 




2 3 4 5 6 

time [s] 



(a) 



(b) 



2 3 4 5 6 

time [s] 



(c) 



Fig. 8. Simulation results — undirected graph: (a) Trajectory of the 8 vehicles: the communication topology is the cycle graph 
Cs (the initial position is marked with a circle and the final position with a diamond; a star indicates the switching time); 
(b) Time evolution of the x-, j/-coordinates of the vehicles (top and bottom, respectively): the same color convention as in (a) 
is adopted here; (c) Time history of parameter s. 



5. 1 Undirected communication graph 



following switching signal (see Fig. 8(c)): 



Consider a team of n single-integrator agents, 

t>t{t) = Vt{t), i e {l,...,n}, 

where p,(t) = [Ptx{t) , p^y{t)]^ G R^ and v>,{t) e R^ 
denote respectively the position and the input of vehicle i 
at time t. Let the control input of agent i be of the form, 

u,it) = (s^ - l)p,(t) + s Y, (P.;-W - SP.W), (12) 

where A/'(i) denotes the set of vertices adjacent to vertex 
i in the communication graptPl. Then, the collective 
dynamics of the group of vehicles adopting control (12), 
can be written as, 

p(t) = (-A(s) ® l2)p(t), 

where p = [pf , . . . , p^]"^ e E,^" and "0" denotes the 
Kronecker product. 

Fig. 8(a) shows the trajectory of 8 vehicles implementing 
the control law (12), when the communication topology 
is the cycle graph Cg (the vehicles are initially on the 
vertices of a regular octagon centered at the origin, and 
their position is marked with a circle in the figure) . For 
the sake of illustration, in our simulation we selected the 



m 



Note that we implicitly assume here that parameter s is 
broadcast to all the agents by a supervisor, via a centralized 
transmitter. 



for t e [0, 1) sec, 
-1 for t e [1, 8] sec. 



The time evolution of the x-, y-coordinates of the agents 
is reported in Fig. 8(b). As it is evident in Figs. 8(a) 
and 8(b), the vehicles first converge towards the origin by 
maintaining equal interdistances, and then the even and 
odd agents cluster in two distinct groups (recall Prop. 2). 

5.2 Directed communication graph 

A scenario similar to that described in Sect. 5.1 is con- 
sidered in Fig. 9. In this case, 5 single-integrator agents 
implement the control law (12) and communicate using 
a directed graph. This leads to an overall closed-loop 
system of the form, 

p{t) = i-A{V{s))^h)p{t). 

The results in Figs. 9(a)-(c), are relative to the di- 
graph in Fig. 7(a), i.e., D5. Fig. 9(a) shows the tra- 
jectory of the 5 vehicles starting from the position 
p(0) = [-1.5, -3, 2.5, 0, -4, 2.5, 4, -3, 1.5, 3]^. 
As shown in Fig. 9(d), 



s{t) 



l/cos(167r/5) for t e [0, 50) sec. 



1 



for t G [50, 100] sec. 



Fig. 9(b) displays the time evolution of the x-, y- 
coordinates of the agents. From Figs. 9(a) and 9(b), we 
can see that the vehicles first move counterclockwise 
along a common elliptical trajectory with frequency 
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Fig. 9. Simulation results — directed graph: (a) Trajectory of the 5 vehicles: the communication topology is the digraph in 
Fig. 7(a), (the initial position is marked with a circle and the final position with a diamond; a star indicates the switching 
time); (b) Time evolution of the x-, j/-coordinates of the vehicles (top and bottom, respectively); (c) Time history of param- 
eter s; (d) Trajectory of the 5 vehicles: the communication topology is the digraph in Fig. 7(b); (e) Time evolution of the 
X-, y-coordinates of the vehicles; (f) Time history of parameter s. 



/(5) = ^tan(167r/5) and phases 0,(5) = ^^^j^j^, 
i e {1, . . . , 5}, and then rendezvous at the point, 

5 ..5 

(-E?'-(50), -^p,,(50)), 

i—l ?'— 1 

i.e., they achieve average consensus (recall Prop. 14). 

Finally, the results in Figs. 9(d)-(f), are relative to the 
digraph in Fig. 7(b). Fig. 9(d) shows the trajectory of the 
agents with p(0) = [-4, -3, 2.5, 0, -4.5, 1, 4, -3.25, -2, 3]^. 
In this case (see Fig. 9(f)), parameter s evolves 
according to: 



m = 



-1.6889 for t e [0, 50) sec, 
1 for t e [50, 100] sec. 



From Figs. 9(d) and 9(e), we observe that the vehi- 
cles first move counterclockwise on five closed orbits, 



and then rendezvous at a single point (recall item 2 in 
Example 2). 

It is worth observing here that the two simple instances 
described in this subsection, look promising building 
blocks for more sophisticated multi-agent tasks, such as, 
containment control (Ji et al., 2008) and cooperative pa- 
trolling (Pasqualetti et ai, 2012). 

6 Conclusions and future work 

In this paper we have presented a generalization of the 
standard consensus protocol, called deformed consensus 
protocol, and wc have analyzed its stability properties in 
terms of the real parameter s for some special families of 
undirected and directed graphs. Preliminary results for 
arbitrary graph topologies are also provided: however, 
some work still needs to be done in order to precisely 
characterize in graph-theoretical terms, the variegated 
behavior of the deformed consensus protocol (a glimpse 
of such a richness of behaviors is provided by Exam- 
ples 1 and 2). The proposed theory has been illustrated 
via extensive numerical simulations and examples. 
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In future works, we aim at studying the properties of the 
deformed consensus protocol when the (weighted) com- 
munication graph is not fixed but changes over time, and 
at investigating other "parametric" Laplacian matrices 
besides the deformed Laplacian (a first step towards this 
direction has been recently done in (Morbidi, 2013, sub- 
mitted)). 
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